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O ' Abstract 



RG flow of central charge Ceg is investigated for the two boundary sine-Gordon model 
at the free Fermi limit. Thermodynamic Bethe ansatz approach is used to check the 
non-monotonic decreasing properties of Ceg, its resonance, and the modification of c^e due 
Qh! to the mismatch of the periodicity of a Lagrangian parameter and that of the boundary 
scattering parameter. The detailed analysis uses the singularity structure of the system 
on the complex rapidity plane. 
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H ■ 1 Introduction 



The c-theorem states that effective central charge Ceg decreases monotonically as the 
bulk system size increases [T]. This c-theorem is best studied and confirmed in the periodic 
bulk system of integrable field theories (IFT). IFT provides a non-perturbative study of 
the scale dependence of Cefr through thermodynamic Bethe ansats (TBA) [21 El- TBA is a 
set of integral equations for spectral densities of the system, which uses the kernel derived 
from the scattering amplitude. Ceg is written in terms of the weighted integration of 
particle energy. It is remarkable that IR properties obtained from scattering matrices can 
be used to find the RG-fiows from UV to IR. There are many examples which demonstrate 
the success of the TBA approach [H El El 13 E| • 

Notwithstanding the success, Ceg does not obey c-theorem always. In the presence 
of open boundaries non-monotonic behavior has been expected as well as resonances jH]. 
In addition, in some cases TBA does not fully reflect the UV properties of Lagrangian 
Pm ITTl 1121 Uni- This originates from the singularity structure of the system. 

It is to be noted that even though the TBA equations are defined over the real rapidity 
variables, the physical rapidity space is not restricted to the real rapidity line. For a 
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periodic bulk system the rapidity is extended to complex values with — vr < Im(6') < tt. 
The singlarity structure on the complex rapidity space provides all the information of 
particle spectra including bound states. 

If the system has a boundary, then the physical rapidity plane reduces to — 7r/2 < 
Im^^ < 7r/2 and RG flow of c^s can be seen in terms of R-channel TBA Depending 
on the boundary parameter ranges, the RG flow of Ces is known to be much affected by 
the existence of singularities near the real rapidity axis. Especially, TBA may not have a 
finite solution at a certain parameter value, where singularities sit on the real rapidity. In 
this case, TBA is not well-defined beyond the parameter range and should be modified so 
that UV properties from the Lagrangian matches with IR properties from the scattering 
matrix. This demonstrates that complete understanding of the singularity structure is 
essential for understanding RG flow of Ceg. 

In this paper, we study the singularity structure of the boundary sine-Gordon model 
(bsG) on a strip. bsG is given as |9j 

" / "^'^ / {^(^'^V^)'-2/i(cos(2M-l)| 

/OO POO 
dy cos(60L - Xl) - / dy cos(6</)i?, - xr) ■ 

-oo J —oo 

where (f>L{y) = V^(0,?/) ( (t>R{y) = '^{R^y)) is the boundary field living at the left (right) 
edge. R is the width of the strip representing the size of the system. The coupling constant 
6^ is real and is restricted to be less than 1. (&^ is scaled by Svr from the conventional 
choice = Svrfe^). 

In fact, the integrability allows additional boundary terms, 

r^y d^dy) am H^ydMy) 



TT dy TT J_^ dy 

Unlike in the one-boundary action, one cannot eliminate the whole term by introducing 
the bulk term [T^ . 



ah r°° 
/ dx dyd^dy(p{x,y) 

Jo J-oo 



- 1 dy'^'^^^y^ 



Therefore, the general action of bsG on a strip will be 

pR POO f -[ ^ ni h 

A = dx -(9,y.)2-2Mcos(2M-l) -— / , 

Jo J-oo I4vr J TT dy 

/OO poo 
dy cos(60L - x^''^) - 2fif / dy cos(60^ - x^""^) . (1-1) 
-oo J-oo 

We will concentrate on the free fermionic limit of bsG, = 1/2 (or = 47r) ^ni- This 
limit is chosen for simplicity and clarity since, even though (R-channel) TBA is known 
when the bulk scattering is diagonal, there appears a complicated singularity structure 
and imposes instability in the numerical analysis even in the massles case [T^ I13j. 
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The purpose of this paper is to understand RG-flow of Ces and its dependence on 
the Lagrangian parameters through the singularity structure of the system. In section |21 
we briefly summarize the relevant TBA analysis and analyze the domain of parameters. 
According to the boundary pole structure, the boundary parameters are classified into 
four Regimes: First, Lagrangian parameter x is divided into two domains: < x < 
and its extended one 7r/2 < % < vr. Then, each domain is divided according to the 
scattering parameter 77 appearing boundary scattering matrix in ()2-5|) below: < rj < tt /2 
and 7r/2<r7<7r. In section 01 massless limit of the theory is considered and the 
boundary parameter dependence of Ces is summarized. Massless limit reduces the number 
of independent parameters and simplify the analysis a little bit. In section |3J massive 
theory with small boundary strength at both edges (called NN-type) is analyzed. In 
sectional massive theory with large boundary strength fis at both edges (called DD-type) 
is analyzed. In sectional the case with small /x^ at one edge and large /i^ the other edge 
(called ND-type) is analyzed. Section [7| is the conclusion. In Appendix Ising model is 
briefly presented. 

2 Massive theory and boundary conditions 

The TBA of the free Fermi limit is simply given as 

eo = 2rcosh^ (2-1) 

with M = 27Tfi and r = MR. This does not mean that the system is trivial. All 
the essential difficulties lies in the non-trivial boundary action, which induces complex 
singularity structures. The effective central charge is given as 

Ceff = — / de cosh e log Z{e) =c- 24A , (2-2) 

J -00 

where 

z{e) = 1 + Xo{e) e-^« + Xdie) e-^'° . (2-3) 

Here Aq and are fugacity due to the boundary effect, compactly written as jTTj 

Ao(^) = WKt" + KtKt + WKn" + K^K-^' = TrjK^KR) 
XM = {Kt^Kr - Kt~Kl+){K++K],- - K+-K^+) = Det(Kii^«) , 

where the overlines denote for complex conjugate operation. Subscript L and R stand 
for the left and right edge contribution. Kl,r is obtained from the boundary scattering 
amplitude Rl,r{9), K{0) = R{i^ — 6). Explicitly they are written in terms of soliton and 
anti-soliton basis [Qj, 

K^^{n,^\e) = p^{r],^\e) G{r],^\e) , K^'-^{7],^\e) = Q^{v,^\e) G{r],^\e) (2-4) 

where ^ = ± and 

P^{i], 'd\6) = —i sinh(^^) cos{r]) cosh('(9) — ^ cosh(6') sin(r7) sinh(79) 



3 



le 



-^2ia ginh(^) cosh(^) 



(2cos(2±^)cos(2^)) (2cosh(2^)cosh(2^)) ' 
After some manipulation, and Aq are written as 



(2-5) 



A„ 



tan tan tann tann 



2 / V 2 

^ /Vr - ^ fVR + 1-^ 
tan tan 



tanh ( — ^ tanh ^ ^ ~^ 



(2-6) 



Ao = 2 ^ sinh^ 9 cos rji cos r/^ cosh -(9^ cosh + cosh^ 9 sin rj^ sin rjji sinh -(9^ sinh 

+ cosh^ 9 sinh^ 9 cos(2a)) Gi(r/, ^, u) Gr{ti^ ^, u) . (2-7) 



(2-8) 



From the fugacity one can check the periodic property of Z in 

Z{9 + i27T) = Z{9) , Z{9 + in) = Z{9) e^'^^^^ / \q{9) . 



K^~^ in (j2-4p represents the sohton number violation scattering at the edge and 
contains the paremeter a which is identified with the Lagrangian parameter a in (jl-lj) . 
This is because the boundary derivative term in shifts the conjugate momentum of 
the field 11 = 0/(27r) by —a5{x — R)/tt, and as the consequence, the phase of the soliton 
operator in |17; is shifted by a, e*°. 

The scattering parameters t] and d are related to the Lagrangian parameters rriB = 
/iBV^27r/M and < X < TT [TF: 

cos(?7/2) cosh('i9/2) =mBCOsx, sin(?7/2) sinh('(9/2) = mBsm.x- (2-9) 

7] has the mirror symmetry with respect to x = ^/2, v{x) = ~ X)i whose explicit 
relation is plotted in Figured 




Figured f] v.s. x 
Note that the boundary bound state spectrum is 120] 

Sb = M cos u, u = 7] — 7^/2 or tt — [i] — tt/2) 
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(2-10) 



where < u < 7r/2. Thus, if < r/ < 7r/2 there is no boundary bound state. According to 
the existence of a boundary bound state, we conveniently classify the boundary parameters 
into four Regimes: 



Case < X < 7r/2 : when 7t/2 < rj < tt Regime I 

when < rj < 71 /2 Regime II 



Case 7r/2 < X < TT : when 7t/2 < rj < n Regime III 

when < Tj < 11/2 Regime IV. 



It is convenient to define Xc which satisfies ri{xc) = vr/2. Xc is restricted to the value 
< Xc < 7r/2. When < < 1/V2 (which will be called N-type boundary condition), 
Xc = and only Regimes I and III are allowed. On the other hand, when 1/V2 < 
(called D-type boundary condition), < Xc < 7r/4 and all of the four Regimes are allowed. 

One may wonder if Regime I and III (II and IV) are the same each other. However, this 
is not the case. The reason is that scattering data and Lagrangian data do not coincide 
each other in Regime I and III (II and IV) : The periodicity of ?7l/_r in the scattering 
amplitude differs from that in the Lagrangian. rj^/pt in ()2-5|) is periodic in vr and the 
symmetry x —X restricts the domain of x as 

< X < vr/2 in TEA. (2-11) 

On the other hand, x in Lagrangian ()1-1|) is periodic in 27r and is restricted to 

< X < ^ in Lagrangian. (2-12) 

The periodicity of x from the scattering data does not match with the one given from the 
Lagrangian. This suggests that Ces obtained from TEA approach is in trouble when 
X > 7r/2 (Regime III and IV) since TEA uses the IR scattering data. 

One would simply double the period of x in (12- 5|) . Eut this does not solve the problem. 
Suppose one extends the zone of x of one edge to < x < ti"- Then, Cefr(x) = Cefr(7r — x)- 
On the other hand, replacing x into vr — x in the action will change the sign of the 
boundary term since fiB will be replaced with —fis- This means that c^s is insensitive 
to the sign of the boundary action and will have a cusp at x = 7r/2, which is not sound 
from the physics point of view TT^. Thus Ces given in ()3-3|) has to be to modified in the 
extended parameter space. This possibilities are to be treated in detail in next sections. 

The parameter a has no periodicity problem: a in ()2-5|) is periodic in vr. On the 
other hand, a-term in the action (jl-l|) will give an exponential factor in the partition 

function e"^*"^—^. Since 6A0o/(27r) is an integer (winding number), the periodicity of 
a is vr, same as in the scattering data. Thus, the range of a may be reduced to the domain 
—71/2 < a < 71 /2 . Using the symmetry a —a in Aq, the domain can be further reduced 
to 

0<a<7r/2. (2-13) 
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3 Massless limit 



The effect of the boundary parameters on Ces is easier to understand if the number 
of parameters can be reduced. The massless theory is a good example of this. In this 
section, we summarize the singularity structure at the massless limit of bSG. 

The massless limit is obtained by rescaling the parameters approriately ^Hl- Intro- 
ducing a large rapidity parameter ^o, one rescales the mass scale {M/2)e^° —>■ M, rapidity 
6 ^ 6 + 6q, and boundary parameter ^ + Oq. Therefore, one has the parameter 
relation, 

e'^ = 47r/i|/M, r] = 2x. (3-1) 
The rescaled quantities simplify the scattering data of ()2-4|l : 

^ ^ 2 cosh (2^) ^ ' 2 cosh (2^) ^ ' 

The effective central charge is written as 

12r 

Ce« = — dee'logZ{e) (3-3) 

^ J —00 

where Z{9) is the massless limit of ()2-3|) with eo = 2re^ and 

cosrje^ 2 + cos(2a) e 



Ao 



2cosh(2^)cosh(''«-^^ 



A. = tanh (^) tlnh (^) . (3-4) 

In this massless limit, only the combination of Xl and Xr enters in c^s with 1] = 2x = 
2{xr — Xl) and boundary bound state disappears in ()3-2|) . This simplifies the analysis 
greatly. 

The fugacity is written in a simple form when both of the boundaries are given in an 
extreme condition so that they can be treated as a conformal field theory (CFT). When 
both edges are in Dirichlet limit {'^l = "^r ~^ 00) or in Neumann limit {'^l = "^r ^ —00): 
Ao = 2 cos ( and A^ = 1 with ( = 2x for Dirichlet and ( = 2a for Neumann condition. In 
this limit 

Ceff = 1 - SC/n^ , (3-5) 
where — vr < ( < tt. This result is obtained using the integration result for 



POO 

/ (islog(l + 2cosCe^ 
Jo 



' 6 2 



The conformal dimension of the ground state is given as A = CV(87r^)> which reflects 
the duality between the role of a of Neumann condition and the role of x of Dirichlet 
condition in CFT of a compact boson theory. 

When one edge is given in Neumann limit and the other edge in Dirichlet limit, Aq = 
and \d = —1. Then the effective central charge is independent of x and a, 

Ceff = -^, (3-6) 
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since 



/•oo 2 

The conformal dimension of the ground state is given as A = 1/16. 

The effect of a on Ceg is plotted in figure |21 at x = 0. The dependence of Ceg on a 
is manifest when both boundaries are given as Neumann condition; the smaller i) is, the 
greater the effect. 




_2l 1 U 

pi/4 pi/2 



Figure 121 Ces v.s. a at xl = XR = 0. 

The figure El implies that when plotting Ces v.s. with fixed a and x there will appear 
a bump at a sufficently large a. The bump is the result of the interference between the 
ground state and an excited state. This will be discussed in detail in section jH 
On the other hand, the branch-cut singularities arise in c^s when 

Z{e) = . (3-7) 

This type of singularity will be called Zo-singularity {ZqS). There are infinite number of 
ZoS"s and all of them lie at Im(6') = ±7r/2. When rj = tt = vr/2) and a = 0, however, 
there appears ZqS at the real rapidity axis since 

ZW=(l-ta„h(^i^))(l-ta„h(^?£^)). 

This is the source of trouble for the mismatch of the period of x in fwo approaches, TBA 
and Lagrangian. One cannot simply extend the domain ofx<7'"/2tox>^/2 but has 
to take care of the singularity at x = 7r/2 properly jTUl IHl 1121 US] • 

One way to proceed is to analytically continue x < by complexifying rj and rotating 
around at vr so that one can arrive at other domain of x > 7r/2. Under this operation, 
one notices that for the ground state a singularity at m/2 + Op [Op is real and positive 
crosses the real rapidity line and sits at —m/2 + 9p (or vice versa) fT^. This singularity 
will be called Zo-crossing singularity {ZqCS). 
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ZqCS pushes the integration contour around the cross singularity position so that Ceg 
is put as 

deeUogZ{e). (3- 



Ceff 



24r n 12r 
e^- + ^ 



TT 



The first term in RHS is due to the singularity crossing and the second one is the real 
axis contribution where represents the principal value of the integration. 

Figure IHl is the plot obtained from ()3-3|) when x < ^/2, and ()3-8j) when % > 7r/2. In 
this modified version, the curve is not cusped at x = ^/2, but is smooth and the period 
of CefT matches with the period of Lagrangian. 




a=0, i3=50 
(x=ji/2, i3=0 
ci=ji/4, i3=0 
-X- a=0, ^0 



Figure 13 Ces v.s. x for < x < '^r from the modified Ces- 



4 NN-type boundary condition 

In this section, the boundary strength itlb oi both edges are considered to be small, 
^^,-R) ^ \j\p2 (called NN-type). In this case the boundary bound state exists and 
Regime I and III (7r/2 < 77 < vr) are allowed. 
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In figure El Ces is plotted against log(r) for Regime I wlien a = 0. In UV limit 
(r 0), Ceff —>■ 1 and is independent of x- This limit corresponds to the Neumann limit 
of the massless theory. In IR limit {r ^ 1), Ceg — >■ 0. In between UV and IR, Cefr is not 
monotonically decreasing but there appears a resonance bump. This resonant bump is 
enhanced as m"^^^^ approaches 0, which is seen in figure El 

The resonance implies that the ground state interferes with the excited state, the 
boundary bound state. This is explained from the singularity structure. Note that a 
boundary bound state provides a pair of poles in Z ()2-3|) . which will be called Z^o singu- 
larity (^oo'S'), 

Z = oo. (4-1) 
ZoqS lies at ±i(7r — rj) with < (tt — r^) < 7r/2. This is shown in figure (HI 




-8 -6 -4 -2 2 4 6 8 

Re(9) 



Figure El Contour plot of Z^o singularities. 
Crossing singularities lie at ttie imaginary 
axis in the physical strip. 




Re(8) 



Figure 13 Contour plot of Zq singularities. 
Crossing singularities are squeezed 
between Z^o crossing singularities. 



Most of ZoS"s lie at ±(z7r/2 -^Oo)-, where > is real and depends on the scale as 
well as the boundary parameters. At Oo = 0, there is a ubiquitous singularity pair. In 
addition to this, an additional pair exists at ±i6i\f in this NN-type boundary condition 
where < 9^ < {tc — rj). As r 1, 9^ ^ {n — rj). The singularities are given in figurejTl 

In this NN-type, the ground state with itlb can be degenerate with the ground state 
with —niB as r — i> oo and — 0. In the intermediate scale, the two states are not 
completeley degenerate but interfere with each other. This is because as x ~^ or 
niB — 0, ?7 — s> vr and Z^CS at ±i(7r — rf) moves close to the real rapidity line. At the 
same time ZqCS sqeezed between this Z^oCS moves close to the real hne too. In fact, 
ZqoS at ±z(7r — rj) and ZqS at ±.i9N are crossing and the ground state with — /i^ can be 
obtained from the ground state with /xb by taking care of the crossing singularity. (See 
Eq. ()4-2|) below). As the singularities lie near the real rapidity axis the two states interfere 
each other strongly. This resonance effect is first noticed in Ising system in the presence 
of weak magnetic field in p. (See Appendix for Ising model). 

In Regime III (x > 7r/2), one has to count in crossing singularities, which cure Ceg 
to have the right periodicity of x- When x = 7r/2, = vr and the singularity sits at 
^ = 0. Thus X = 7r/2 is the singular point in the parametric space. One may analytically 
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continue x iiito the extended domain 7r/2 < x < vr. (Equivalently, one may analytically 
continue tub around = 0. See this possibility in boundary Lee- Yang case [TU] ) 

Under the analytic continuation, Z^CS as well as ZqCS is crossing at the same 
time. The ground state contribution is obtained by minimizing the number of crossing 
singularities and taking the crossing singularities close to the imaginary axis. 

In figure ini and m the the crossing singularities are pointed with arrows. Z^CS is 
given as ±i(7r — rjfj) and ZqCS as ±i9j\f. When only one edge's x exceeds 7r/2, Cefr is 
modified as 

Ceff = -— (sin^^v - sinr^^j + — j dOe' \ogZ{e) , (4-2) 

where the right edge is chosen to have xr > ^/2- On is the real function of r, m^, x, and 
a. Even though the exact dependence of 6'jv is not known analytically, it can be calculated 
numerically. When r ^ oo, 9 ^ ^ rjji. When r ^ 0, ^ 9o and 6*0 vanishes as — > 0. 
At UV limit, Ces = 1 and at IR limit, Ces = 0. In figure |H1 plotted is Ces v.s. log(r) when 
^^/^) _ g_2^ anfi a = xl = 0. 
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-6 -5 -4 -3-2-1 1 2 

log(r) 



FigureEl CeS v.s. log(r) at = = 0.2 and a = xl = 0. Curves with xr = 0, 7r/4 
belong to Regime I, and XR = 37r/4, vr to Regime III. 



In figure |H1 rn^B ~ used. When ^ one needs to calculate Ceg in two 

ways and chooses the minimum as the ground state result among the two : One with 
(Xl, Xr) and the other with (xl, xr) where Xl = vr - ^XL^and xr = - Xr- This is 
because bsG with {xl^Xr) is equivalent to one with {xliXr) under the field operation, 
bip ^ bip + 71 and —<p), but the crossing singularity contribution can be different. 

Suppose both x's exceed 7i/2, (e.g., xl > 7r/2 and xr > 7r/2). Then one may rearrange 
both parameters as x < 7r/2 so that the crossing singularities do not contribute. If one 
would count in the crossing singularities of both edges, then Ces would be of an excited 
state result rather than of the ground state one. 
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The parameter a also affects Ccs- Especially, at UV limit the effect is manifest and 
Ceif 7^ 1 when a 7^ 0. One may convince oneself that at UV limit Ceg v.s. a reproduces 
the CFT limit ()3-5|) . This behavior is plotted in figure IHl (Refer to figure IT^ below also). 




-6 -5 -4 -3 -2-1 1 

log(r) 



Figure ini Ces v.s. log(r) at XR = XL = and niB = 0.001 with various a. 

At IR limit, however, a effect vanishes and c^s —>■ 0. Therefore, RG-flow between UV 
and IR is not monotonically decreasing. A resonance bump is also seen at the intermediate 
scale due to the singularities near the real axis. 

5 DD-type boundary condition 

In this section, both edges are restricted to DD-type (mj^^'^^ > 1/^2). In DD-type, all 
the four Regimes are to be considered. To simplify the analysis, we put rriB = m^B^ = m^^\ 




Figure [TUl Cefr v.s. log(r): 



m 



(R) 



1000 and XL = a = 0. 
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At IR limit Ccs ^ 0. At UV limit Ccs —>■ 1 when a = 0. However, at the intermediate 
scale, the plot Ces v.s. r (or Ces v.s. log(r)) shows a new feature. In figureHUl Ces v.s. log(r) 
is plotted for = 1000 and a = 0. (Here xr = 0, vr/S belong to Regime II, and xr = 
3tt/S, 7i/2 to Regime I.) There appears a plateau around log(r) ^ —8 (here ~ 15 3> 1), 
which corresponds to the massless limit : Ceff v.s.x at this plateau reproduces the value 
in the massless limit figure El The plateau appears around — log(r/2) ^ i)/2 ^ 1, which 
is between — log(r/2) ^ and — log(r/2) ^ 1. The bigger tub, the wider the plateau 
becomes. This is why the massless limit is easily obtained when ^ 1. The Dirichlet 
limit of the massless theory corresponds to the extreme limit rriB oo. 

Note that the massless limit does not coincide with the UV limit. For a sufficiently 
large ^ 1, provides a large characteristic energy scale. One may rescale the rapidity 
around to see the massless limit, which is considered in Sec. El UV limit appears when 
the energy scale is greater than (length scale is smaller than — log(r/2) ^ 

In regime I (xc < X < ^/2), the boundary bound state exists. Thus Z^o has the 
singularity in the physical strip at ±i{n — r^L//?). In figure ITTl Z^qS is shown when a = 0, 
rriB = 1000, Xl = 0, = 27r/5 and log(r) ^ —10. (Here rjji = 47r/5 and rji = 0). 




ZqS lies at lm{d) = ±7r/2 as r <C 1. The singularities accummulate around at ^ = 
±(i7r/2 ± log(r/2)) along with the isolated one at ^ = ibz7r/2 as shown in figure IT2l 
However, as r ^ 1, a pair of (crossing) singularities appear at the pure imaginary rapidity 
axis ±i9N with 9^ real, {n — rj) < 9^ < 7r/2. 

Regime HI (7r/2 < x < (tt — Xc)) is the analytically continuation of Regime I and Ces is 
modified according to the crossing singularities. Suppose Xl = and n/2 < xr < (tt — Xc)- 
We need to consider two cases seperately, r large and r small. When ZqCS is given as 
ib(^2? =t ^7r/2) (the case when r is small), 

24r / \ Qr f°° 

CefT = -— cosh ^D- sin r/fij + — y d9 e'^ log Z {9) . (5-1) 

When Zq crossing singularity is given as ii^Ar (the case when r is large), 

24r / \ 6r /""^ 

Ceff(x) = sin^^ - sinr/^) + — 4 d9 e' \ogZ{9) . (5-2) 
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In UV limit, ()5-H) gives Ces 1- In IR limit, a.s 9^ 
behavior is seen in figure IT^ when xr = 37r/5. 



VR (El, Ceff ^ 0. This 




Figure El Ces v.s. logr (left) and Cefr v.s. r (right) with 
in Regime III and XR = 47r/5, vr are in Regime IV. 



4" 



1000. Xr = 3vr/5 is 



In Regime II, (0 < x < (^r — Xc)), there is no boundary bound states and therefore 
no ZooS on the physical strip. ZqS is given as Im(6') = 7t/2. Analytically continuing 
to Regime IV ((tt — Xc) < X < t'"), one has ZqCS at ±{9c, ± i7r/2) with 6'^) real and 
positive. It is noted however, even though there is no boundary bound states, Z^oCS is 
at ±i(7r — Tjn)-, which is outside of physical strip. Thus, Ceg is given as the same form of 
(IS3D: 



24r 



CefT 



vr V 



( COsh^n — 



D - Sm TjR 



6r 



rf^e^ logZ(^). 



(5-3) 



The r-dependence is given in figure IT^ for xr = 47r/5 and tt. In the intermediated 
scale, the plateau reproduces the massless hmit in figure El As r 1, one sees that Ceg 
has the tendency of decreasing linearly in r In fact, at IR limit and as a result 



24r 

Ceff(x) = 

TT 



l_sinr^^)+0(e-2^) . 



(5-4) 



The linearity of Cefr in r indicates that if one changes the sign of one edge, the free energy 
density /(r) in IR limit becomes bigger than the case with the same sign: 



A/ 



TT 



2AR 



Aceff = M - M smr]R. 



(5-5) 



The additional free energy corresponds to the one soliton mass M compensated by the 
right edge boundary excitation energy M smrj^. Thus, one may subtract this IR amount 
from Ceff and equivalently, re-normalize the free energy density. 

Finally, a effect is vanishingly small at IR limit. However, the a effect is enhanced at 
r smaller than the plateau region. Especially, at UV limit Cefr is not 1 but depends on 
values of a. Cefr v.s. a at UV limit reproduces the CFT limit, which is plotted in figure 
d 
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Figure ITU Ceg v.s. a: r = exp(— 23) and ms = 100. 



6 ND-type boundary condition 

In this section, we will give a few comments when combination of N-type at one edge 
and D-type at the other edge is given. Suppose the boundary strength at the left edge is 
neglectingly small (m^'' ~ 0) and that of the right edge is large (m^"* ^ l/y/2). In this 
case the singularity crossing effect can be neglected. (Note that one may choose x < 7r/2 
if = using the action symmetry). Still, there remains a plateau, correponding to 
the massless limit in addition to the UV limit. 

At UV, one sees the enhanced role of the compact boson property, c^s ^ 1 — 12(a/7r)^, 
independent of x- At the plateau, one may confirm that c^s is independent of x and a 
and approaches to —1/2. This results in A = 1/16, the dimension of the ground state. 
One may see this behavior in figure IT^ 





Figure El Cgfr v.s. log(r) with various a: Figure El c^s v.s. log(r) with various x- 



m 



(L/R) 
B 



(10-2, 10^) and xl/r = (0, vr/4) . m^^ ' = (10"^ 10^) and a = 



At IR Ccs —>■ 0. In addition, as r S> 1, there apppears a resonance bump. The 
resonance behavior is more enhanced as xr ^/2, since the singularities move close to 
the real rapidity axis. This is clearly seen in figure El 
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If the boundary strength ttib is not in the extreme case, the singularity crossing con- 
tribution is not significant. When one edge isinO < x < '^/'^ ^ind the other edge in 
71 / 2 < X < ^! using the action symmetry one may re- arrange the N-type boundary in 
7r/2 < X < TT (Regime III) and D-type boundary in0<x<7r/2 (Regime II) so that one 
counts in the singularity crossing from the N-type boundary, which will give the ground 
state contribution. 

7 Conclusion 

To summarize, we consider the free fermion limit of the massive boundary sine Gorodon 
model on a strip and investigate the boundary parameter effect on the RG flow of Ces- For 
a certain range of parameters, the flow shows Ceff = at IR limit, and Ces ^ 1 — 12(Q;/7r)^ 
at UV limit. However, depending on the paramter ranges the RG fiow is not monotonic 
and and shows much more diverse effects at the intermediate scale. All the anomalous 
behavior originates from the singularity structure on the complex rapidity space. 

In NN-type boundary condition {m^^^^ < l-\/2), UV hmit and massless limit coincide. 
In the intermediate scale a resonance appears due to the singularity close to the real 
rapidity axis. If one edge isinO < x < ^/'^ (Regime I), and the other edge is in 
7r/2 < X < TT (Regime III), then the singularity crossing is to be counted to modify the 
Ceff formula. 

In DD-type boundary condition [m'^^^^ > l-s/2), UV limit does not coincide with 
the massless limit. The massless limit appears as the plateau in c^s v.s. log(r), which 
is characterized by the scale r ~ e~'^ . The Dirichlet limit of this plateau reproduces the 
massless result, Ceg 1 — 12(x/vr)^ with X = Xr~ Xr- 

If one edge isin0<x<7r/2 (Regime I or II) and the other edge isin7r/2<x<7r 
(Regime III or IV), then the singularity crossing modifies the Ces formula. Especially 
when one edge is in Regime II (0 < x < 7r/2 and < 77 < 7r/2), and the other edge in 
Regime IV (7r/2 < x < tt and < 77 < 7r/2), then is linearly decreasing in the scale r, 
which indicates the opposite sign of boundary terms excites the one-particle state of the 
same sign of boundary terms. Similar behavior is obtained for Ising model when opposite 
boundary magnetic fields are applied to both edges. The magnetic field should not be too 
strong (—1 < k < 0), otherwise the opposite fields spoil the field theory and result in the 
complex free energy. 

In ND-type, (m^^ < Iv^ and m^"* > 1V2), still UV limit and massless limit have 
their own domain. In the extreme limit {m^^^ and m^"* ^ 1), one may see that 
UV limit is the compact bosonic CFT with Ceg ^ 1 — 12(Q;/7r)^. On the other hand, the 
massless limit gives Ces — > —1/2. 

The similar singularity effect on RG flow is expected to be persistent when bsG goes 
beyond free Fermi limit. However, the presence of more species of particles and boundary 
bound states makes RG flow more complicated and the numerical analysis will be more 
difficult due to the nemerical instability. 

It is noted that the singularity analysis for the case of opposite boundary terms will be 
useful to study other boundary fleld theories. In Ising model, the magnitude of opposite 
boundary flelds is limited since the fleld theory breaks down. This does not happen 
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in bsG. On the other hand, in boundary affine Toda field theory, the boundary actions 
are given in 3 kinds only. Among them the (-)boundary action is not well understood 
[221 123 1211 • seems that the complete understanding of the boundary action requires the 
proper study of the singularity structure at the presence of (-)boundary action condition. 
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Appendix 

Ising model's reflection amplitude is given as 

Rkid) = ita.n\i{— - -] — — 7-1 

where k = 1 — h"^ with h the boundary magnetic fleld. (Note that /i^ is scaled by 2M 
where M the particle mass so that h is a dimensionless number.) Thus 

:r^/^N . , A; + cosh6' ,„ 

and the fugacity becomes 



2 ki + cosh^ kn + cosh^ 



X = KJ<, = Unk'[-)^—^^l—^^. (7-3) 



The effective central charge is given as 

Cefr = — / d9 cosh e log Zi3i„g(^) = c - 24A , (7-4) 

where e = 2rcosh^ with r = RM. Zising(^) = (l + A(6') e""*^^-') satisfies the periodicity, 
^ising(^ + ^27r) = Zising(^) and Z^ngiO + ZTi) = Zising(^) e^°(^) /Aq (^) . 

According to singularity structure, the domain of parameters are classified into three 
regimes: 0</c<l, — 1<A;<1 and k < —1. In each regime, the magnetic field h and 
—h is to be distinguished. In fact, /i = (or A; = 1) is the singular point in the parameter 
space h since one of ZqS lies at real rapidity axis. The domains h > and h < are 
connected through the branch singularity. Thus two domains are not equivalent due to 
the existence of crossing singularities. 

Regime < A; < 1 is the weak magnetic field regime {0 < < 1) similar to the 
N-type in the text. There is a boundary bound state, which makes a pair of Z^oCS at 
±2^ where < ^ < 7r/2 and cos^ = k. A pair of ZqCS at ±i^jv with < < ^ is 



16 



squeezed between the pair of ZooCS. Ces is modified when the boundary magnetic fields 
are of opposite sign {hihfi < 0), 

Ceff = sin 6'Ar - sin n + — 4 d6 e log ^ising(6') • (7-5) 



3.5 




-6 -5-4-3-2-1 1 2 3 4 



log(r) 

Figure [T71 Cgfr v.s. log(r) for < k < 1 (weak boundary magnetic field). 



In regime —l<k<0 (boundary magnetic field 1 < h"^ < 2), there is no boundary 
bound state. Nevertheless, ZooCS exists at ±i(^ + 7r/2) (0 < ^ < vr/2), outside of the 
physical strip. All ZqCS lie at Im(^^) = 7i/2 and ZqCS is given as ±{6d ± ivr/2). When 
the boundary magnetic fields are of opposite sign, 

Ceff = ( cosh 9d - cosH + ^ 4 dOe log Zising(6') . (7-6) 

TT V / TT^ 

In figure UHl Ceff with hihn > is plotted. The case with hihn < is contrasted with 
that of same sign in figure IT^ where the curves with the same sign are nearly overlapped 
and are not distinguished each other. In all cases, there is no resonance. 




Figure ITHl Ceff v.s. log(r) with — 1 < < 0. 
Boundary magnetic field of both edges are 
with same sign. 



- \= hp= sqrt(5)/2 
\^ hp= sqrt(6)/2 

- \^ hp= sqrt(7)/2 

- h|^= -hp= sqrt(5)/2 
_ h|^= -hp= sqrt(6)/2 
_ \^ -hp= sqrt(7)/2 




Figure IT^ Ceff v.s. log(r) with — 1 < < 0. 
Boundary magnetic field of both edges are 
with opposite sign 
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Note that in figure ^1 Ces decreases linearly in r as r ^ 1 and when opposite boundary 
fields are given. (At IR limit, 6d = and the integrated part vanishes in ()7-6p so that 
the linear coefficient is proportional to (1 — cos.^).) This trend of linearity is clearly seen 
in figure QUI 




Figure QUI Ceff v.s. r for — 1 < < 0. 



The linearity of Ces in the Regime — 1 < < shows the similarity to D-type of bsG. 
However, unlike bsG the plateau does not appear in this Regime but in Regime k < —1 
(boundary mag netic field > 2). As given in figure |^ the stronger the boundary field 
becomes, the wider appears the plateau region which is identified as the massless limit. 
The massless limit is distinguished from the UV limit even though c^s 1/2 in both 
limits. When k > —1, the massless limit coincides with the UV limit. 



0. 



0. 



0. 



0. 



0. 




h|_=hp.3O0 
ti|_=hj^.3O00 
^ h|_=hp=30000 



-0.1' I I ' ' ' ' L 

-30 -25 -20 -15 -10 -5 

log(r) 

Figure OTl Cgfr v.s. log(r) for k < —1. 



When k < —1, the opposite boundary magnetic field makes the field theory unstable. 
The breakdown is due to ZooCS at ±(^ ± ivr), where cosh^ = \k\ with real ^ > 0. 
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The singularity lies outside of the physical strip and becomes 'resonant poles'. This 
crossing singularities result in a complex-valued Ces and thus, free energy dendsity becomes 
complex valued : 



Ceff = cosh 9d + i sinh U + ^i d9 eUog Zising(^) , (7-7) 

where ZqCS is at ±{6£, ±i7r/2). 

Finally, if one edge is given with a vanishing boundary magnetic field k ^ 1 and the 
other edge with a strong one k —oo, the plateau appears corresponding to CFT limit 
of the mixed boundary condition, whose Cgfr = — 1 and the ground state has the conformal 
dimension 1/16. (See figure |221) 




^ h|^=0.01, hp-3000 
_0_ h^=0.01, hp-30000 



-30 -25 -20 -15 -10 -5 

liig(r) 

Figure 1221 Cefr v.s. log(r). A weak (strong) boundary field is given at one (the other) edge. 
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